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Abstract 


Interstellar SiO was discovered shortly after CO but it has been detected 
mainly in high density and high temperature regions associated with outflow 
sources. We present a new model of interstellar silicon chemistry that ex- 
plains the lack of SiO detections in cold clouds and contains an exponential 
temperature dependence for the SiO abundance. A key aspect of the model 
is the sensitivity of SiO production by neutral silicon reactions to density and 
temperature, which arises from the dependence of the rate coefficients on the 
population of the excited fine structure levels of the silicon atom. This effect 
was originally pointed out in the context of neutral reactions of carbon and 
oxygen by Graff, who noted that the leading term in neutral atom-molecule 
interactions involves the quadrupole moment of the atom. Similar to the 
case of carbon, the requirement that Si has a quadrupole moment requires 
population of the J = 1 level, which lies 111K above the J = 0 ground 
state and has a critical density > 10 6 cm -3 . The SiO abundance then 
vhas'a temperature dependence proportional to exp(-lll/ T) and a quadratic 
density dependence for n <n„ . As part of the explanation of the the lack of 
SiO detections at low temperatures and densities, our model also emphasizes 
the small efficiencies of the production routes and the correspondingly long 
times needed to reach equilibrium. Measurements of the abundance of SiO, 
in conjunction with theory, can provide information on the physical proper- 
ties of interstellar clouds such as the abundances of oxygen bearing molecules 
and the depletion of interstellar silicon. 

Subject headings: interstellar: abundances - interstellar: molecules - molecu- 
lar processes 



I. INTRODUCTION 


Interstellar SiO was detected at mm wavelengths (Wilson et al. 1971) 
shortly after the discovery of CO (Wilson, Jefferts, and Penzias 1970). Since 
then, CO has been observed to be the most widely distributed molecule 
other than H 2 in the interstellar medium. In contrast, SiO has been detected 
mainly in regions with high density and temperature associated with outflows 
from young stellar objects. Although silicon is likely to be strongly depleted 
on dust grains in dense clouds and SiO requires a higher density to excite 
than CO, the almost complete lack of detections in quiescent regions cannot 
be understood solely in these terms. In this paper, we explain the lack of 
extensive SiO detections in terms of the sensitivity of its production by neu- 
tral silicon reactions to density and temperature, an effect directly associated 
with the excitation of the atomic fine structure levels (Graff 1989). 

Interstellar silicon chemistry was first discussed in detail by Turner and 
Dalgarno (1977); other early work was carried out by Solomon and Klemperer 
(1971), Prasad and Huntress (1980) and Millar (1980). Silicon chemistry in 
shocks has been studied by Hartquist, Oppenheimer, and Dalgarno (1980) 
and more recently by Neufeld and Dalgarno (1989). Circumstellar silicon 
chemistry has been investigated by Scalo and Slavsky (1980) and Clegg et 
al. (1983). Recent theoretical and experimental studies of silicon reactions 
and new astronomical observations now provide the stimulation to develop 
further the interstellar chemistry of silicon. 

From radio observations it is apparent that the interstellar chemistry of 
silicon differs from that of carbon, oxygen, nitrogen, and sulfur. The ele- 
mental abundance of Si is twice that of S and only an order of magnitude 
less than C, N, or 0. Yet these species occur in many molecules in widely 
distributed interstellar molecular clouds whereas gas phase silicon is detected 
only as SiO and SiS (Morris et al. 1975, Irvine et al. 1986) and then only in 
hot dense regions perturbed by star formation. (These two molecules are also 
observed in the circumstellar envelopes of highly evolved red giants; SiC and 
SiC 2 have recently been detected in IRC +10216 (Thaddeus, Cummins, and 
Linke 1984; Cernicharo et al. 1989)). To address the issue of the dependence 
of SiO on temperature, Ziurys, Friberg, and Irvine (1989) searched for its 
J = 2 — 1 emission from both cold as well as hot molecular gas. Despite 
the fact that the observations of the colder regions were made down to very 
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low noise levels, SiO was detected only in sources with kinetic temperatures 
greater than 30 K. Because the dipole moment of SiO (3.1 Debye) is much 
larger than that of CO, the SiO emission arises from dense regions, typically 
with n(Hj) > 10 4 — 10 5 cm -3 . Moreover, in the sample considered by Ziurys, 
Friberg, and Irvine (1989), all but one of the SiO detections were made in 
or near outflows from young stellar objects and in the supernova remnant 
IC 443. To determine the abundance variation with temperature, Ziurys, 
Friberg, and Irvine (1989) compared the SiO column densities with those of 
HCN, another molecule requiring high density excitation and found primar- 
ily in dense regions but not manifesting any strong temperature sensitivity. 
Over the temperature range from roughly 15 K to 200 K, they found that the 
relative abundance of SiO varies roughly as exp(-90/ T). 

The observations of Ziurys, Friberg, and Irvine (1989) have prompted us 
to find a theoretical basis for a temperature depedent SiO abundance. No 
such temperature sensitivity is known in interstellar carbon chemistry and 
earlier silicon chemical models do not contain any strong temperature effects 
at low temperatures. However, recent theoretical work by Graff (1989) shows 
how the neutral reaction rates of carbon and oxygen atoms, which depend 
on the population, ordering, and spacing of the atomic fine structure levels, 
could have a strong temperature dependence. In the case of carbon, with a J 
= 0, symmetric, ground state, strong interactions with other neutral systems 
require significant population of the excited states, the first of which lies at 
23.6 K above the ground state. The resultant temperature variation of the 
carbon chemistry may not be very pronounced because of the small energy 
level spacing relative to the kinetic temperature of clouds and also because of 
the dominance of many other chemical pathways. Oxygen, with opposite fine 
structure, does not show this behavior. However, silicon has the same fine- 
structure ordering as carbon but with larger separations: the first level is at 
111 K and the second at 321 K. In this paper we propose that SiO in shielded 
regions is formed mainly by neutral reactions of Si atoms in the 7=1 level 
with OH and O 2 . By applying Graff’s ideas, we find that the interstellar 
SiO abundance contains an exponential temperature dependence. We also 
find that the SiO abundance is sensitive to density because the fine structure 
levels of Si are likely to be be sub-thermally populated at densities less than 
10 4 — 10 s cm -3 . Our results explain the absence of SiO detections in cold, 
dark clouds by Ziurys et al. Whether or not the the predicted temperature is 
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the same as reported by these authors is less certain because the conditions 
for the applicability of our model may not be satisfied for the sources they 
discussed. 

II. SILICON CHEMISTRY 

Although much basic laboratory and theoretical information on astro- 
physically interesting silicon compounds is lacking, we are in a better position 
than earlier investigators to evaluate the various reaction pathways. Some 
of the key neutral reactions were actually measured more than a decade ago 
(Husain and Norris 1978, Swearengen, Davies, and Niemcyzk 1978). Im- 
portant recent progress includes the measurement of silicon ion-molecule re- 
action rates (e.g. Bohme, Wlodek, and Fox 1988) and the development of 
techniques for the extrapolation of room temperature rates to lower temper- 
atures (e.g., Clary 1988, Graff 1989). Advances have also been made in the 
thermal-chemical data base. For example, the data for silicon hydrides (Cur- 
tiss and Pople 1988 and references therein) imply that all reactions of Hj with 
both neutral and charged mono-silicon hydrides with up to four H atoms are 
endothermic. Hence the first-stage species in the ion-molecule chemistry of 
silicon, SiH + and SiH, , will not be hydrogenated further and may in fact be 
destroyed by molecular hydrogen. Even if they are not destroyed in this way, 
they provide relatively weak pathways to heavier silicon-bearing molecules, in 
contrast with interstellar carbon chemistry where hydrides provide efficient 
routes for the synthesis of CO and other molecules. 

The most important feature of our chemical model is the population de- 
pendence of the reaction rates of atomic silicon with the neutral species, OH 
and Oj. Graff (1989) has recently investigated neutral atom reactions at very 
low temperatures, which are dominated by the long-range interaction of the 
atomic quadrupole with the dipole or quadrupole moment of the reactant 
radical or molecule. The low temperature behavior of these reaction rate 
coefficients is sensitive to the population of the ground state fine structure 
levels, which determines the atomic quadrupole moment. Graff illustrates 
the effect by comparing reactions of 0 and C atoms with the radicals OH 
and CH. Because these atoms have opposite fine-structure orderings (Table 
1), their low-temperature reaction rate coefficients are very different: oxy- 
gen will interact strongly with neutral radicals and molecules and carbon 
weakly below 20 K. Graff suggests that the same principles apply to neu- 
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actions of other atoms important in the interstellar medium and we list some 
of the relevant physical parameters in Table 1. We omit nitrogen because 
its (2p) 3 4 Si/ 2 ground state hats no fine structure. Because its quadrupole 
moment also vanishes, we speculate that N reacts the weakest of all the 
abundant heavy atoms unless other dynamical effects govern its low temper- 
ature reactivity. 

We extend Graff’s ideas to the silicon atom reactions with radicals and 
molecules that are important in interstellar clouds. The fine structure level 
ordering of silicon is similar to carbon, but the spacing is larger: E\q = 111 K 
and E n — 321 K (Table 1). Thus we expect these reaction rates to be 
proportional to exp(-lll/!T) at low temperature. Carbon atom reactions 
will be governed by a similar population factor but, because of the smaller 
level spacing, significant effects can only be expected for cold clouds with 
T < 15 K. The silicon atom reactions also depend on density because the 
fine structure level population in interstellar clouds may not be in thermal 
equilibrium. For the temperatures of interest here, T < 200 K, the effects 
of the J — 2 level can be ignored and statistical equilibrium theory gives 
the following result for the the total steady state fractional population of the 
first excited level of Si 

Pi = 3 (n/[n(l + 3 exp(-Tu>/T) + n^]) exp(-T 10 /T) (1) 

where Tier = j4io/cio is the critical density for the J = 1 — 0 transition, 
assumed to be optically thin, is the spontaneous decay probability, and 
cio is the collisional de-excitation rate coefficient for the J = 1 — 0 transition. 
Recent calculations of the de-excitation rate coefficient of the corresponding 
collisions of He atoms and H 2 molecules with 0 atoms (Monteiro and Flower 
1986, Flower 1989) indicate that cio < 10 -11 cm 3 s -1 for T < 100 K and 
Cio << C 21 , C 2 o- In the absence of similar calculations for Si, we use the 
results for 0 as a guide. If the H 2 + Si rate coefficient ci 0 = 10 -11 cm 3 s _1 for 
T < 100 K, the critical density for the J = 1 — 0 transition is 8 x 10 6 cm -3 . 
The corresponding critical density for the J = 2 — 1 transition would be 
at least 10 times smaller and might give interesting physical effects at low 
densities and T > 100 K. 

Table 2 lists the key reactions in our model of the silicon chemistry of 
interstellar clouds. We assume that the clouds are well shielded from ex- 
ternal stellar and interstellar ultraviolet radiation. We list the adopted rate 
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coefficients and give references to laboratory measurements if they exist. 
Otherwise we follow standard practice and conventional wisdom in estimat- 
ing rate coefficents. If the rate coefficient depends on temperature, we give 
the exponent n if it follows a power law (in T n ) or indicate by “exp” that 
the rate has the exponential temperature dependence exp(-lll/!T) associated 
with silicon atom reactions. For brevity we omit all endothermic reactions. 
They are unimportant for the low temperature regions under consideration 
here but they are crucial for shock and other high temperature chemistry. In 
the following paragraphs, we discuss the background for the several classes 
of reactions listed in Table 2. For ion-molecule reactions with neutrals with 
large dipole moments, we use n = —0.6, following the discussion of Clary 
(1988). 

Neutral Reactions 

We consider three neutral reactions in the synthesis of silicon molecules 
from Si and the most abundant oxygen species other than CO (for which the 
relevant reactions are endothermic). Although the reaction of Si with H 2 0 
is exothermic, we expect it to have a large enough activation energy to make 
it unimportant at low temperatures. The primary neutral reactions for the 
synthesis of SiO are reactions 1 and 2 of Table 2, 

Si + OH -* SiO + H 
Si -f- 0 2 — * SiO 4- 0 


Only the rate coefficient for reaction 2 has been measured and here the two 
room temperature experiments give very different results. Swearengen et 
al. (1978) suggest that the larger value, obtained with a flash photolysis 
technique is incorrect because of heating or the production of other reactive 
species. In this paper we adopt the smaller rate coefficient. Theoretical 
considerations support this choice. Because 0 2 does not have a permanent 
electric dipole moment, the leading term in the Si-0 2 interaction is of the 
quadrupole-quadrupole type, which leads to a reduction in the Langevin rate 
coefficient by an order of magnitude (Graff 1989). Both reactions 1 and 2 will 
be reduced at low temperatures by the effects discussed above and expressed 
by equation (1). 

Applying the results of Graff (1989) to Si, we write the rate coefficients 
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for these reactions as 


fci = ifci(oo)pi (2) 

where ki(oo) (i = 1,2) is the value at high temperatures (relative to the 
fine-structure splitting). We have taken the fraction of the J — 1 states 
that contributes to the reaction to be 1/3 and we have ignored any residual 
interactions associated with the J = 1 — 0 state, again following Graff. Even 
with the smaller rate coefficient, reaction 2 is probably the dominant SiO 
production process because of the large abundance expected for Oj relative 
to OH. There is only one fairly weak route to Si02 (reaction 3, SiO + OH) 
because the reaction of SiO and O 2 is endothermic. 

Cosmic Ray Induced Photodestruction 

Although neutral reactions occupy the central role in our chemical model, 
cosmic rays are also important. In addition to directly ionizing hydrogen and 
helium, they generate UV radiation by inelastic electron scattering from H 2 
which photoionizes Si and dissociates silicon molecules (Prasad and Taraf- 
dar 1983; Sternberg, Dalgarno, and Lepp 1987; Gredel, Lepp, and Dalgarno 
1987). The last group of authors finds that a cosmic ray ionization rate ( 
(defined for H 2 ) leads to destruction rates for interstellar molecules that are 
of order 10 2 — 10 3 £, depending on the photodestruction cross section. We 
have made estimates for silicon species using the U V emission spectrum given 
by Gredel et al., which spans the range from 900- 1650 A. We have estimated 
the rate for Si, with ionization threshold at 1520 A, by scaling the result of 
Gredel et al. for C, using the Si photoionization cross section (Hudson, 1971) 
and taking the grain albedo as 0.5. The result, Gc U » = 3 x 10 -14 s -1 is 3 times 
larger than for C because the Si photoionization cross section matches better 
the strong UV emission by H 2 near 1500 and 1250 A. If the albedo is taken 
to be 0.9, the rate would be 5 times larger. Our estimates for SiO and Si02 
are much less certain because we know little more than the photodissociation 
thresholds of these species, 1510 and 2650 A, respectively. 

Ion-Molecule Reactions 

Pathways to SiO formation can be initiated by radiative association of 
Si + with H 2 (reaction 8) and by the two exothermic reactions of Si + with 
OH and H 2 0, (reactions 9 and 10 in Table 2). It is noteworthy that all 
of the measured ion-molecule reaction rates listed in Table 2 are an order 
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of magnitude less than the Langevin rate. Despite this fact, we ignore any 
possible activation energies as is customary in dealing with the low tempera- 
ture behavior of ion-molecule reaction rates. Only an upper limit is reported 
for reaction 9 and it can be ignored. The thermal-chemical data needed to 
determine whether the reaction SiOjJ" -f H2 — ► Si02H + + H (reaction 12) is 
exothermic is missing at this time. The corresponding reactions for CO^ and 
SOj do occur (Anicich and Huntress 1986). However, this point is not very 
important here because the production of SiO^" is probably weak. 

The most important ions for destroying SiO and Si02 axe C + and He + . 
We have introduced characteristic branching ratios for the various channels 
that arise when He + destroys SiO and Si02- Many reactions of He + with 
other neutral molecules have been measured (Anicich and Huntress 1986) 
and serve as the basis for the following estimates. Because charge exchange 
is generally found to be weak (and frequently does not occur with oxides), 
we argue that these branchings are small (B 3 and 64 < 0.1). Measurements 
for CO and NO, suggest that, for SiO, Si + rather than 0 + is favored (B 2 ~ 1 
and B\ < 0.1). Similarly, data for CO2 and SO2 suggest that the dominant 
branch for Si02 is SiO + but that Si + and 0 + may be produced 20 % of the 
time. For numerical estimates, we adopt 61 = 0.8, b 2 = 0.1, b 3 = 0.05, and 
64 = 0.05 and Bi = 0.8, B 2 = 0.2, and B 3 = 0.1. 

Charge Exchange and Neutralization 

Charge exchange with neutral atoms is potentially an important destruc- 
tion mechanism for Si + but there is considerable uncertainty in the rate 
coefficients at the temperatures of interstellar clouds (e.g., see the review by 
Dalgarno 1988). We have followed the guidelines and estimates of Pequignot 
and Aldrovandi (1986) in assuming that Si + charge exchanges rapidly with 
the following heavy atoms: Mg, Fe, Al, Ca, and Ni. Of course all of these 
elements are likely to be strongly depleted in grains, thereby reducing the 
effect. However, Si + can be efficiently neutralized by the population of nega- 
tively charged solid particles and large molecules such as polycyclic aromatic 
hydrocarbons expected to exist in shielded interstellar clouds (e.g., see the 
review by Omont 1986). Without commiting ourselves to the exact identity 
of this component of the interstellar medium, we have included their neu- 
tralization of Si + as reaction 28 in Table 2, denoting them by the symbol 
PM~. The relatively small abundance of such “particles” is offset by the 
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large cross section for neutralization, and this process is likely to dominate 
charge exchange with neutral atoms. 

in. CHEMICAL MODEL 

The production pathways for the chemistry discussed in § II are repre- 
sented schematically in Figure 1. The dominant pathways for the synthesis 
of SiO and Si0 2 involve the reactions of Si with OH and 0 2 and Si + with 
H 2 0. In this section, we derive approximate, analytic formulae which exhibit 
the dependence of the various silicon species abundances on temperature, 
density, and cosmic ray ionization rate. 

As mentioned in § II, the production of silicon hydrides is extremely in- 
efficient. The SiH + produced by proton transfer from H3 to Si does not 
react with H 2 and most of it recombines to Si again. The only other signifi- 
cant route to SiH, radiative association (reaction 8) followed by dissociative 
recombination of SiHj, is also weak, as we will show below. Therefore we 
initially neglect the silicon hydrides and their ions because they are relatively 
unimportant as progenitors of SiO and Si0 2 , except at very low temperatures. 

The SiO + ion has negligible abundance because it reacts rapidly with H 2 
(reaction 11) to form SiOH + . It is unlikely that SiOH + interacts further with 
H 2 ; instead it recombines to form SiO. Thus the interaction of Si + with OH 
leads to the following ionic pathways to SiO: 

Si + + OH - SiO + + H 
SiO + + H 2 SiOH + + H 
SiOH + + e -> SiO + H 

The corresponding neutral pathways are initiated by reactions 1 and 2. 

In the approximation of ignoring the silicon hydrides, we solve for the 
steady state silicon abundances using the three balance equations for Si + , 
SiO, and Si0 2 and the conservation condition for all gas phase silicon. The 
steady state equations for the fractional abundances of species X, ®(X) = 
n(X)/n(H 2 ), formed from species Y and Z are written in the form 

Dx(X) = Px( Y)x(Z) 

where D and P are appropriate destruction and production factors. The 
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balance equations are given below and the P’s and D’s are listed in Table 3: 


£>ix( Si + ) = P 4 x(Si) + P 5 x(SiO) + P 6 x(Si0 2 ) (3) 

D 2 x(SiO) = Pix(Si) + P 2 *(Si + ) + Pjx(Si0 2 ) (4) 

£> 3 x(Si0 2 ) = Prx(SiO) (5) 

Solving these equations yields the following ratios among the silicon species: 

x(Si0 2 )/x(SiO) = P 7 /D 3 (6) 

x(SiO)/x(Si) = [Pi + P 2 P 4 /D 1 ]/ A (7) 

where 

A = [D 2 - P 2 P 5 /Pi] - [P 7 /D 3 }[P 3 + P 2 P 6 /i?i] (8) 


IV. RESULTS 

We use the formulae derived in § III to estimate the abundances of the 
silicon species as a function of the relevant astrophysical parameters. The 
essential nature of the temperature and density dependence can be most 
clearly understood by first assuming that Si0 2 and SiH can be neglected, 
which we later show to be a good approximation. The factor A in the 
denominator of equation (8) then reduces to A = D 2 — P 2 P$/Di and the 
second term can be neglected as long as x(H 2 0) < 10~ 4 . Consequently the 
SiO/Si abundance ratio can be written in a form where the contributions to 
molecule formation by neutral reactions with Si and ion reactions with Si + 
occur additively: 

x(SiO)/x(Si) = R f exp(— 111/T) + Z (9) 

where 

R = Pi(oo)/Z> 2 (10) 

Z = P 2 P 4 /D 1 D 2 < \ Pi/Di (11) 

/ = n/[n(l + 3exp(-lll/T)) + n^j. (12) 

The factors R, Z, and / depend only weakly on temperature. The symbol 
Pi ( 00 ) is the production factor Pi in Table 3 evaluated at 300 K. Taking into 
account the conservation condition for silicon yields 

x(SiO)/x Si = [R f exp(— 111/T) + Z}/{ 1 + Rfexp(-lU/T) + Z] (13) 
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where xsi is the total abundance of all gas phase silicon species. This equation 
has an exponential temperature dependence for T\ < T < T 2 if Rf is not 
very large and 2 < 1, We can estimate 7\ and T 2 by ignoring the term 
3exp(— 111/T) in equation 12 and obtain 

T\ = 111 f[ln(Rf/Z)] (14) 

T 2 = lU/[ln(Rf/(l + Z))\ (15) 

if Rf/( 1 + Z) > 1; otherwise T 2 = 00 . Substituting typical values of the pa- 
rameters into R, f , and Z, expressing them in convenient units, and assuming 
that Cio = 10 -11 cm 3 s -1 , these quantities become 

R ~ 16 n B [ 35 ( 02 )— s + 0 . 1 x(OH)_ 8 ]/£_i 7 (16) 

/ ~ ns/(ns -(- 8.2) (17) 

Z ~ 0.0025 x(H 2 0)_ 6 /[x(PM-)_ 7 r-° 5 + 0.25x(A)_« + 0.005 x(H 2 0)_«] 

(18) 

where we have introduced the notation, n 5 = »/ 10 s cm 3 , etc. 

According to equation 13 and the approximate forms for R , Z, and /, the 
abundance of SiO is sensitive to the abundances of OH, O 2 and H 2 0. At this 
time there is little observational information concerning the abundances of 
O 2 , and H 2 O because of the extreme difficulty in observing the emission of 
these molecules from the ground. Indeed, the apportioning of oxygen among 
the various oxygen-bearing species in molecular clouds is widely considered to 
be one of the outstanding problems in interstellar chemistry. Radio observa- 
tions of OH (e.g., Crutcher 1979, Irvine et al. al. 1986) indicate a range from 
10 -8 — 10 -7 for the OH abundance in clouds with densities n(H 2 ) < 10 4 cm -3 . 
Theory suggests that the OH abundance decreases with density (Langer and 
Graedel 1989) and we adopt an intermediate value for estimation purposes, 
x(OH) = 3 x 10 -8 . Emission from molecular oxygen has not yet been de- 
tected and the smallest upper limit, x(0 2 ) < 5 x 10~®, has been obtained 
with the isotopic variant, 16 0 18 0, for the p Oph cloud (Liszt and Vanden 
Bout 1985, Goldsmith et al. 1985). Here we adopt a slightly smaller value 
z(02) = 3 x 10~ 6 representative of the results of chemical models. Radio ob- 
servations of H 2 O (and its isotopic variants) involve transitions from highly 
excited states in hot, dense regions (n( H 2 ) > 10 4 cm -3 , T > 10 2 K) and sug- 
gest that x(H 2 0) ~ 10 -6 . A possible detection of infrared absorption from 
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the ground state of H 2 O leads to x(H 2 0) < 3 x 10 -7 — 2 x 10 -6 for the 
material along the line of sight towards Orion-BN (Knacke, Larson, and Noll 
1989). Because H 2 0 can be expected to freeze out on grains at relatively high 
temperatures ~ 100 K, we expect the H 2 0 abundance to vary from region 
to region and, for estimation purposes, assume that it is in the range from 

10- 8 - 10- 5 . 

The quantities Rf and Z are the relative fractions of SiO formed by 
neutral and ionic reaction pathways, respectively. One or the other must 
approach unity in order for a substantial fraction of silicon to be converted 
into molecules in steady state. It is important to note that Rf and Z have 
different dependences on density: Rf ~ n 2 for n < and Z independent 
of n. In order to have efficient production of SiO at high temperatures by 
neutral reactions, the density must be > 5 x 10 4 cm -3 for the choices made 
above for the abundances of oxygen- bearing molecules. The quantity Z is a 
sensitive function of the H 2 0 abundance and of the abundances of the species 
that neutralize Si + . Even for a:(H 2 0) as large as 10 _B , Z is only ~ 0.01 for 
T = 100 K. Because Z is generally small, SiO is produced efficiently by 
neutral reactions in warm, dense clouds, but the role of neutral reactions is 
significantly reduced in cool clouds by the Boltzmann factor, exp(— 111/T). 
More specifically, equation (14) gives the condition that the temperature 
must satisfy in order for neutral reactions to dominate ionic reactions, i.e., 
T > Ti, where T\ is determined approximately from equation 14 using the 
quantity 

RSjZ _ 6-7 x [T- x(PM”)_t 

For example, if x(H 2 0) ~ *(0 2 ), then I\ ~ 25 K. Because H 2 0 is expected 
to freeze out on grains at low temperatures, it is reasonable to expect that 
neutral reactions will continue to dominate production down to very low 
temperatures. Thus, if a:(H 2 0) ~ 0.1x(O 2 ), T\ is reduced to 15 K. 

In Figure 2 we plot the solution to equation 13 for the ratio x(SiO)/*si 
as a function of temperature (plotted as 1/T) for several values of Rf and 
Z. As expected from the above discussion, the fractional SiO abundance 
decreases approximately as exp(— 111/T) if Rf is not too large and Z is 
small. In the case where Rf is large, the exponential decrease does not occur 
until low temperatures, e.g., it starts at T = 50 K for Rf = 10. For Z ^ 0, 
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the exponential decrease disappears at low temperatures and the solution 
goes over to a constant. In the case Rf = 1 and Z = 10 -3 , this changeover 
occurs below 20 K. However, we show below that time dependent effects 
become important at low temperatures and this leveling out probably does 
not occur. 

Although the contribution of SiH to SiO production is usually small, it 
can become important at low temperatures because SiH can interact strongly 
with ground state oxygen atoms. The primary path for SiH formation is 
radiative association of Si + with H 2 (reaction 8) followed by recombination 
with electrons (reaction 24). It is destroyed primarily by 0 (reaction 4) as 
long as x(O) > 10 -6 . Under these conditions the SiH abundance is given 
approximately by the following expression, 

x(SiH) ~ 5 x 10- 3 *(Si + )/x(O)_ 4 . 

Because they are both proportional to the abundance of Si + , we can easily 
compare the efficiency of SiO production by SiH with that of the ionic path- 
way involving Si + and H 2 O (reaction 10). Production by SiH needs to be 
considered only if x(H 2 0) < 2 x 10 -7 . 

It is important to examine the time required to synthesize the silicon 
molecules. It is apparent from the above solutions that efficient conversion 
of silicon atoms into molecules only takes place at high densities. Thus little 
SiO should be produced in the early, low density stages of cloud evolution 
when n < 10 4 cm -3 . This conclusion can be brought into clearer focus by 
considering the time for neutral reactions to produce SiO, which can be 
approximated by 


tchcm - pi(oo)x(0 2 ) + Moo)*(OH)) n f exp( — 111/T)] 1 


and can be written in convenient units as 


tchem = 3 x 10 4 exp(lll/T) 


(1 + 8.2/715) 

n 5 [x(0 2 )-s + 0.1 x(OH)- s 


yr 


For T > 100 K, the time required to achieve equilibrium is fairly long 
at low and intermediate densities, e.g., approximately 3 x 10 7 yr for n = 
10 4 cm~ 3 and x(02) = 3 x 10 -6 and x(OH) = 3 x 10 -8 and 3 x 10 5 yr if 
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the density is increased to n = 10 5 cm -3 . Densities greater than 10® cm -3 
are required before the neutral reactions can produce SiO in a time short 
compared to cloud evolutionary times but, even in this case, the clouds must 
be warm. The exponential increase of the chemical timescale suggests that, 
at low temperatures, where ions dominate its production, the abundance 
of SiO will be very small. Consequently, we believe that it is unlikely that 
much SiO is produced at low temperatures because, not only are the chemical 
pathways inefficient, but the time to produce near-equilibrium abundances 
is much longer than the lifetimes of the dense, cold regions of interstellar 
clouds. 

The effects of non-equilibrium chemistry are illustrated in Figure 3 where 
we show approximate time dependent solutions for x(SiO)/xsi a function 
of temperature (plotted as 1/T) for Rf = 1 and for n — 5 x 10 4 cm -3 , 
®(0 2 ) = 3 x 10 -6 , and x(OH) = 3 x 10 -8 . We give solutions at times 
ranging from 10 6 to 10 9 yr as well as the steady state solution for Z = 0. 
For t < 10 7 yr, the abundances are considerably depressed below the steady 
state solution at low temperatures (T < 25 K). For the case Z = 10 -3 
(ionic chemistry included), we also present results at 10 7 years and in steady 
state. Even though this steady state solution departs from an exponential 
(and becomes constant) below 20 - 25 K, the time dependent solution remains 
nearly exponential down to about 15 K, although with a steeper slope than 
the function ezp(— 111/T). In the limit of T — > 0, the ionic chemistry time 
scale becomes t c hem = 3 x 10 8 yr, independent of density, and probably 
longer than any relevant dynamical time scale. Figure 4 present additional 
illustrations of the theory for higher density regions without ionic reactions 
( Z = 0). For n = 10 6 cm -3 (Fig. 4b), large SiO abundances are achieved 
for temperatures down to 15 K within ~ 10 6 years. For times less than 10 6 
years, the abundance has an exponential dependence on temperature. 

Finally, we can estimate the Si0 2 abundance in warm interstellar clouds 
of the type considered in our theory. From the equations in § III we can write 
in convenient units 

x(Si0 2 )/z(SiO) ~ 0.17 n$ x(0H)_8/C-i7* 

Significant amounts of Si0 2 can be produced in warm regions where SiO 
is abundant. Using the OH abundance adopted above, x(OH) ~ 3 x 10 -8 , 
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we obtain Si02/Si0 ratios of order 0.5 n 5 . Although SiC >2 is not directly 
observable by radio techniques, it might be possible to detect its protonated 
ion, HSiOj", which should have a large dipole moment. We estimate that 

HSiOj /Si0 2 ~ 10- 3 . 

V. SUMMARY 

We have developed a model of the gas phase silicon chemistry for shielded, 
quiescent interstellar clouds in which the SiO abundance has a strong temper- 
ature dependence of the type suggested by Ziurys et al. In this theory, neutral 
reactions of Si with O 2 and, perhaps, OH are mainly responsible for produc- 
ing SiO. These reactions occur only for Si in excited fine structure levels and 
require relatively warm conditions, e.g., T > 30 K. Fairly high densities are 
also necessary to significantly populate the first excited fine structure level 
of Si, (n(H 2 ) > 10 6 cm -3 ), in order to overcome its rapid spontaneous decay. 
Ion-molecule reactions are relatively unimportant in shielded regions in pro- 
ducing SiO and, when they occur, proceed with such low efficiency that very 
long times are needed to achieve observationally significant values. Accord- 
ing to our theory, the production rate of SiO at low temperatures and low 
densities is so small that steady state abundances are unlikely to be achieved 
within cloud evolutionary times. 

In our model, the fraction of gas phase silicon in SiO is sensitive to the 
density, the O 2 and OH abundances, the total gas phase abundance of silicon, 
and the critical density for populating the J = 1 level of Si (see equations 13 
and 16-18). Consequently, in conjuction with observational determinations of 
the SiO abundance, it should be possible to use our theory to determine, or at 
least constrain, presently poorly known physical parameters of appropriate 
interstellar clouds such as the depletion of silicon and the abundances of 
O 2 and H 2 O. For example, measurements of the SiO abundance in warm, 
shielded regions could be used to determine the silicon depletion. In the high 
T limit 

and, because Rf must be at least of order unity to produce SiO efficiently, 
the SiO abundance itself measures (to within a factor of a few) the gas phase 
abundance of silicon. Unfortunately, it is not yet possible to make a quanti- 
tative estimate along these lines because the observations discussed by Ziurys 
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et al for warm regions refer to outflow sources or their environs. The abun- 
dance of SiO in these sources is likely to have a significant contribution from 
the reaction of Si + with H 2 0 (equation 10). In this situation, an estimate 
of the SiO abundance can be made in the context of a photodissociation 
region model (Tielens and Hollenbach et al, 1985) by assuming that SiO is 
produced by reaction 10 and destroyed by UV radiation at a rate G(SiO): 

*(SiO) = kio n x(Si + ) x(H 2 0)/(?(Si0) 

The SiO abundance now depends on the depletion factor for silicon, the 
strength of the radiation field, the density, the H 2 0 abundance, and the 
photodestruction cross section of SiO. Entering reasonable values for these 
quantities leads to abundances of the same order given by Ziurys et al. for 
warm, outflow sources, although with considerable uncertainty because of 
the dependence on so many unknown parameters. 

Obtaining constraints for cold, dark clouds is more difficult because the 
chemistry is almost certainly time-dependent. The steady state expression 
for the SiO abundance (equation 13) reduces to 

x(SiO) ~ Z xsi 

in the low temperature limit where ions contribute. Referring to equation 
18 for Z, we see the potential for obtaining a limit on the H 2 0 abundance 
in cold, dark clouds using the upper limits to the SiO abundance estimated 
by Ziurys et al A low water abundance for cold, dark clouds is consistent 
with the expectation that H 2 0 readily freezes out onto grains during the 
course of cloud evolution at low temperatures. However, the nondetection of 
SiO in these clouds could also arise from severe depletion of silicon species 
onto grains or the long chemical time scale for the formation of SiO or a 
combination of all these effects. 

In summary, our theory explains the failure of Ziurys et al to detect 
SiO in dark clouds in terms of the suppression of both the neutral and ionic 
chemistries, the former because of the low temperature and density and the 
latter by the freezing out of H 2 0. The direct confirmation of the exponential 
temperature dependence associated with the excitation of the fine structure 
levels requires further observational studies of warm, shielded clouds. Despite 
the sensitivity of the theory to many presently unknown reaction rate coeffi- 
cients and to the physical parameters and chemical abundances of interstellar 
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clouds, it does demonstrate the importance of neutral reactions for interstel- 
lar silicon chemistry. Our work also illustrates the significance of Graff’s ideas 
on the low temperature behavior of the rate coefficients for the reaction of 
atoms with molecules and radicals. In focusing in this paper on the most 
abundant silicon molecule, SiO, we have ignored several interesting problems 
in interstellar chemistry such as the synthesis of SiS and the abundance of 
the silicon carbides. We have also not attempted to treat the coupling be- 
tween the silicon in the gas and solid phases other than to recognize that 
most of the silicon is incorporated in dust in dense clouds with depletion 
factor varying from cloud to cloud. These simplifications are justified by the 
current, preliminary nature of both the theoretical and observational aspects 
of interstellar silicon chemistry. 

Our theory suggest several directions for future research. In order to 
strengthen the basis of the theory, it would be extremely valuable to remea- 
sure the fundamental neutral reaction, Si + 0 2 — >■ SiO + 0, as a function of 
temperature at and below 300 K. Calculations of the collisional de-excitation 
rate of the Si fine structure levels by H 2 would help quantify the density 
dependence of the chemistry. Measurements of the 0 2 and H 2 0 abundances 
would permit interesting tests of the theory based on the close correlation 
between SiO and the 0 2 and H 2 0 abundances. More sensitive ground based 
observations of different types of interstellar clouds are required to establish 
the density and temperature dependence of the SiO abundance. One of the 
goals of such observations should be to confirm the prediction of this paper 
that the SiO abundance in dense, shielded clouds varies with temperature 
and density, primarily through the population of the excited fine structure 
levels. 

One of us (A.E.G) would like to acknowledge partial support of this re- 
search by a grant from NASA. The authors are grateful for helpful comm- 
ments from Dennis Downes, Paul Goldsmith, Michel Guelin, Chris McKee, 
and Frank Shu. 


16 



Table 1 


Fine Structure Energy Levels and Transition Probabilities® 


Atom 

Eo/k 

Ei/k 

E 2 /k 

Aio 


C 

0 

23.6 

62.6 

7.93 x 10“ 8 

2.68 x lO" 7 

0 

325.9 

228.1 

0 


8.95 x lO" 5 1.70 x lO" 5 

Si 

0 

111.0 

321.1 

8.25 x 10- 6 

4.20 x lO" 5 

s 

825.4 

571.0 

0 


1.40 x lO' 3 3.02 x 10“ 4 


° The units are K for energies and s 1 for transition probabilities. 
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Table 2 


Silicon Chemical Reactions 


Reaction 


Rate Constant 0 T-Dependence 


1. Si + OH — ► SiO + H 

10-i° 

exp 

2. Si 0 2 — * SiO -)■ 0 

9 x 10~ 12 b 

exp 


2.7 x 10- 1Oc 

exp 

3. SiO + OH — ► Si0 2 + H 

10- 10 


4. SiH + 0 -» SiO + H 

lO" 10 


5. CR UV + Si -*■ Si+ + e 

3 x 10~ 14d 


6. CR UV + SiO Si + 0 

3 x 10~ 14<i 


7. CR UV + Si0 2 -> SiO + 0 

3 x 10~ 14 J 


8. Si+ + H 2 -► SiH^ + hu 

10~ 16 


9. Si + + OH — ► SiO+ + H 

< 10~ n e 


10. Si+ + H 2 0 -► SiOH+ + H 

2.3 x 10~ 1Oe 

-0.6 

11. SiO + + H 2 SiOH+ + H 

3 x lO' 10 ' 


12. SiOj + H 2 -> Si0 2 H+ + H 

2 x 10 -1 ° 


13. Et + Si -> SiH+ + H 2 

lO" 9 


14. HCO+ + Si SiH+ + CO 

lO" 9 
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15. H t + SiO -+ SiOH + + H 2 

lO " 9 

- 0.6 

16. HCO+ + SiO SiOH+ + CO 

io - 9 

- 0.6 

17. He+ + SiO -* 0+ + Si + He 

Bi x IO " 9 

- 0.6 

He+ + SiO -► 0 + Si+ + He 

B 2 x IO ' 9 

- 0.6 

He+ + SiO -► SiO+ + He 

B 3 x IO " 9 

- 0.6 

18. He+ + Si0 2 -► SiO+ + 0 + He 

fei x 10~ 9 

- 0.6 

He + + Si0 2 — * Si + -(- 0 2 + He 

b 2 x IO ' 9 

- 0.6 

He + -(- Si0 2 — ♦ 0+ 4 " SiO -(- He 

b 3 x IO ' 9 

- 0.6 

He+ + Si0 2 -► SiO^ + He 

64 x IO " 9 

- 0.6 

19. C+ + SiO -h. CO + Si+ 

10“ 9 

1 

O 

b> 

20 . C+ + Si0 2 -> CO + SiO+ 

10- 9 

- 0.6 

21 . Si"*” -(- e — * Si hi v 

8.7 x IO -12 ® 

-0.65 

22. SiOH+ + e -+SiO + H 

3.0 x IO " 7 

-0.5 

23. Si0 2 H+ + e -Si0 2 + H 

3.0 x 10- 7 

-0.5 

24. SiH t + e — »SiH + H 

3.0 x IO " 7 

-0.5 

25. H+ + SiO SiO+ + H 

10 -i° 


26. H+ + Si0 2 SiO^ + H 

10- 10 


27. Si + + A -+ A+ + H 

lO-iofc 


28. Si+ + PM” -» Si + PM 

2.4 x 10" 8i 

-0.5 
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° The units for cosmic ray induced UV photodestruction are s -1 ; the 
units for all other rate coefficients are cm 3 s -1 . The numerical values are for 
T = 300 K. If the temperature dependence is a power law, the exponent is 
given in column 3; the symbol exp indicates the temperature dependence is 
exp(— 111/ T). 

b Measured by Swearengen et al (1978). 

c Measured by Husain and Norris et al. (1978). 

d The rates for cosmic ray induced UV photodestruction are based on the 
work of Gredel, Lepp, and Dalgarno (1988). 

e Fahey et al (1981) 

* Anicich and Huntress (1986). 

g Pequignot and Aldrovandi (1986). 

h A stands for atoms that charge exchange rapidly with Si, such as Mg, 
Fe, Al, Ca, and Ni. The abundance of these atoms in dense clouds may be 
of the order of 10 -6 . 

* Omont (1986). 
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Table 3 


Coefficients P» and D* 


Symbol 

Expression 

Order of Magnitude® 

Pi 

k\ x(OH) + ^2 sc( 02 ) 

IO -15 

Pi 

kiQ x(H20) 

2 x 10" 16 

P 3 

61 ^i 8 *(He + ) 

3 x 10 “ 19 

Pi 

Qcruv (Si) 

3 x 10 ‘ 19 

Pi 

B\ fcir x(He + ) 

3 x 10 - 19 

Pi 

b 2 fcis x(He + ) 

3 x lO " 20 

Pr 

k 3 x(OH) 

3 x 10" 18 

Di 

k 27 x(A) + k 23 x(PM ) 
x(H 2 0) + 0.5 fcg 

5 x lO " 16 

D 2 fr 

k 3 aj(OH) + k 17 x(He + ) 
+ 5'cruv(SiO) + fc 19 x(C + ) 

4 x 10‘ 18 

D 3 8 

kis x(He + ) + ycruv(Si0 2 ) 
+ k 2 o x(C + ) 

6 x 10~ 19 


° The order of magnitude estimates are based on the following parameters: 
rc(H 2 ) = 10 5 cm -3 , T = 100K and C = 10 ~ 1 T s" 1 ; x(H 2 0) = 10“ 6 , 
x(0 2 ) = 3 x 10- 6 , x(OH) = 3 x 10 ‘ 8 , x(A) = 10 ~ 6 , and x(PM~) = 10' 7 . 

^9c ruv — G„ uv ln{a 2 ). 
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FIGURE CAPTIONS 


Figure 1. Simplified schematic diagram of interstellar SiO chemistry- 
showing the principle production routes and the destruction of neutrals by 
cosmic ray produced UV radiation (signified by the label i/). Ionic destruc- 
tion has been omitted for simplicity. 

Figure 2. Fraction of silicon in the form of SiO in steady state as a function 
of temperature (plotted as 1/T along the lower axis and T along the upper). 
Three measures of the relative strength of the neutral production, Rf, are 
shown and two of the ion contribution, Z. A nearly exponential decrease 
of SiO with decreasing T down to 15 K requires weak ion (Z < 10 -3 ) and 
moderate neutral ( Rf « 1) production, consistent with conditions expected 
in dense regions of interstellar clouds. 

Figure 3. Fraction of silicon in the form of SiO as a function of tempera- 
ture (plotted as 1/T along the lower axis and T on the upper scale) at different 
times for Rf = 1 and no ion contribution (Z = 0); n = 5 x 10 4 cm -3 , x(C> 2 ) = 
3x 10 -6 , and x(OH) = 3x 10~ 8 . Results including ions ( Z = 10 -3 ) are shown 
for 10 7 yr and steady state. Below 30 K, very long times are needed to ap- 
proach equilibrium and, thus, even with contributions from ions, the solution 
looks like an exponential for T > 15 K over times that are relatively long 
compared with cloud dynamical times. 

Figure 4. Similar to Figure 3 for (a) n = 2 x 10 5 cm -3 (Rf ~ 20) and (b) 
n = 10 6 cm -3 (Rf ~ 200). 
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